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ABSTRACT 

Recent observational results for magnetic fields in molecular clouds reviewed by Crutcher (2012) 
seem to be inconsistent with the predictions of the ambipolar diffusion theory of star formation. 
These include the measured decrease in mass to flux ratio between envelopes and cores, the failure to 
detect any self-gravitating magnetically subcritical clouds, the determination of the flat PDF of the 
total magnetic field strengths implying that there are many clouds with very weak magnetic fields, 
and the observed scaling B cx p^^^ that implies gravitational contraction with weak magnetic fields. 
We consider the problem of magnetic field evolution in turbulent molecular clouds and discuss the 
process of magnetic field diffusion mediated by magnetic reconnection. For this process that we termed 
"reconnection diffusion" we provide a simple physical model and explain that this process is inevitable 
in view of the present day understanding of MHD turbulence. We address the issue of the expected 
magnetization of cores and envelopes in the process of star formation and show that reconnection 
diffusion provides an efficient removal of magnetic flux that depends only on the properties of MHD 
turbulence in the core and the envelope. We show that as the amplitude of turbulence as well as the 
scale of turbulent motions decrease from the envelope to the core of the cloud, the diffusion of the 
magnetic field is faster in the envelope. As a result, the magnetic flux trapped during the collapse in 
the envelope is being released faster than the flux trapped in the core, resulting in much weaker fields 
in envelopes than in cores, as observed. We provide simple semi-analytical model calculations which 
support this conclusion and qualitatively agree with the observational results. Magnetic reconnection 
is also consistent with the lack of subcritical self-gravitating clouds, with the observed flat PDF of 
field strengths, and with the scaling of field strength with density. In addition, we demonstrate that 
reconnection diffusion process can account for the empirical Larson (1981) relations and list a few other 
implications of the reconnection diffusion concept. We argue that magnetic reconnection provides a 
solution to the magnetic fiux problem of star formation that agrees better with observations than the 
long-standing ambipolar diffusion paradigm. Due to the illustrative nature of our simplified model we 
do not seek quantitative agreement, but discuss the complementary nature of our approach to the 3D 
MHD numerical simulations. 

Subject headings: cosmic rays, scattering, MHD, turbulence 



1. INTRODUCTION 

Interstellar media are known to be turbulent and mag- 
netized, and both turbulence and magnetic field are im- 
portant for star formation (see Armstrong et al. 1994, 
Chepurnov & Lazarian 2009, Crutcher 2012). The exist- 
ing star formation paradigm has been developed with the 
concept that near-perfect fiux freezing holds, i.e. that 
the magnetic field is well coupled with ions and elec- 
trons in the media (Alfven 1942). This is assumed in 
magnetically-mediated star formation theory, which was 
founded by the pioneering studies by L. Mestel and L. 
Spitzer (see Mestel & Spitzer 1956, Mestel 1966) and 
brought to a high level of sophistication by other re- 
searchers (see Shu, Adams & Lizano 1987, Mouschovias 
1991, Nakano et al. 2002, Shu et al. 2004, Mouschovias 
et al. 2006). According to the theory, magnetic fields 
slow down and even prevent star formation if the media 
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are sufficiently magnetized. The theory assumes that 
the change of the flux to mass ratio happens due to am- 
bipolar diffusion, i.e. to the drift of neutrals which do not 
feel magnetic fields directly, but only through ion- neutral 
collisions. Naturally, in the presence of gravity, neutrals 
get concentrated towards the center of the gravitational 
potential while magnetic fields resist compression and 
therefore leave the forming protostar (e.g. Mestel 1965). 
The rate of ambipolar diffusion for a cloud in gravita- 
tional equilibrium with the magnetic field depends only 
on the degree of ionization of the media. 

The existing theory makes star formation inefficient for 
magnetically dominated (i.e. subcritical) clouds. The 
low efficiency of star formation corresponds to observa- 
tions (e.g. Zuckerman & Evans 1974), which is usually 
interpreted as a strong argument in support of the above 
scenario. This however does not solve all the problems; 
at the same time, for clouds dominated by gravity, i.e. 
supercritical clouds, this scenario does not work as mag- 
netic fields do not have time to leave the cloud through 
ambipolar diffusion. Therefore for supercritical clouds 
magnetic fields should be dragged into the star, forming 
stars with magnetizations far in excess of the observed 
ones (see GaUi et al. 2006, Johns-Krull 2007). 
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Magnetic fields are important at all stages of star for- 
mation. In many instances the ideas of star formation 
based exclusively on ambipolar diffusion have been chal- 
lenged by observations (Troland & Heiles 1986, Shu et al. 
2006, Crutcher et al. 2009, 2010a, see Crutcher 2012 for 
a review) . While the interpretation of particular observa- 
tions is the subject of scientific debates (see Mouschovias 
& Tassis 2009), it is suggestive that there may be addi- 
tional processes that the classical theory does not take 
into account. The primary suspect is turbulence, which 
is ubiquitous in interstellar media and molecular cloud^ 
(see Larson 1981, Armstrong et al. 1994, Elmegreen & 
Falgarone 1996, Lazarian & Pogosyan 2000, Stanimirovic 
& Lazarian 2001, Heyer & Brunt 2004, Padoan et al. 
2006, 2009, Chepurnov & Lazarian 2010, Burkhart et 
al. 2010, 2012). Turbulence has revolutionized the field 
of star formation (see Vazquez- Semadeni et al. 1995, 
2000, Ballesteros-Paredes et al. 1999, 2007, Elmegreen 
2000, 2002, McKee & Tan 2003, Elmegreen & Scab 2004, 
MacLow & Klessen 2004, McKee & Ostriker 2007), but 
the treatment of the turbulent magnetic fields stayed 
within the fiux freezing paradigm. 

The cornerstone concept of magnetic fiux freezing has 
been challenged recently. On the basis of the model of 
fast magnetic reconnection in Lazarian & Vishniac (1999, 
henceforth LV99), Lazarian (2005) claimed that the re- 
moval of magnetic fields from turbulent plasma can hap- 
pen during star formation due to magnetic reconnection 
rather than slow ambipolar drift (see also Lazarian & 
Vishniac 2009). The process that was later termed re- 
connection diffusion does not depend on the degree of 
ionization, but rather on the properties of turbulence. 
The numerical confirmation of the idea was presented 
for molecular clouds and circumstellar accretion disks in 
Santos-Lima et al. (2010, 2012). The numerical test- 
ing of the LV99 reconnection model that is at the core 
of the concept of reconnection diffusion was successfully 
tested in Kowal (2009, 2012). Compared to this test- 
ing of reconnection diffusion, the testing in Kowal et al. 
(2009, 2012) was performed with much higher numerical 
resolution of the reconnection region and much better 
control of turbulence and other input parameters. In ad- 
dition, more recent formal mathematical studies aimed 
at understanding of magnetic field dynamics in turbu- 
lent fluids supported the LV99 quantitative conclusions 
(Eyink 2011, Eyink, Lazarian & Vishniac 2011). More- 
over, the state of observations of magnetic fields in re- 
gions of star formation has progressed significantly (e.g. 
Crutcher 2012), so that comparisons of observational re- 
sults with the theory are possible. These developments 
motivate us to further study the implications of the re- 
connection diffusion concept for star formation. 

We note that the interstellar medium is collisional (Ya- 
mada et al. 2006), and therefore ideas of collisionless 
reconnection (see Shay & Drake 1998, Shay et al. 1998, 
Bhattacharjee et al. 2005, Cassaket al. 2006) are not ap- 
plicable. Even if they were, they would not increase the 
rate of reconnection diffusion, which is controlled by the 
eddy diffusivity. Similarly, the ideas related to the plas- 

The presence of turbulence is an observational fact, while its 
sources are debatable. Supernova explosions, different instabilities, 
e.g. magnetorotational instability, have been discussed in the liter- 
ature. Flows induced by the gravitational instabiltiy (see Vazqucz- 
Semadeni et al. 2011) are also expected to be turbulent. 




Fig. 1. — Upper panel. Sweet-Parker reconnection. A is limited 
by resistivity and is small. Lower panel: reconnection according 
to the LV99 model. A is determined by turbulent field wandering 
and can be large. From Lazarian, Vishniac & Cho (2004). 

modial instability (see Shibata & Tanuma 2001, Lour- 
reiro et al. 2007, Uzdensky et al. 2010, Huang et al. 
2011) do not change the reconnection rates either^ 

In what follows, in §2 we discuss the reconnection of 
magnetic fields in turbulent fiuids, in §3 we discuss re- 
cent observational results reviewed by Crutcher (2012) 
that can be interpreted as contradicting the model of star 
formation mediated by ambipolar diffusion. We provide 
in §4 a very simple semi-analytical model of magnetic 
field diffusion out of magnetized clouds and discuss its 
predictions. In §5 we discuss other observational tests. 
Astrophysical implications of reconnection diffusion are 
briefiy discussed in §6, while the discussion and summary 
are provided in §7 and §8 respectively. 

2. RECONNECTION DIFFUSION OF MAGNETIC FIELD 

First of all, we should discuss how magnetic fields can 
diffuse in a conducting medium in the presence of turbu- 
lence. For this purpose we should consider the process 
of magnetic reconnection. In the absence of magnetic 
reconnection frozen-in magnetic field lines cannot cross 
each other, and the turbulent fiuid should become in- 
creasingly entangled and eventually get locked. Using an 
analogy one can think of the fiuid acquiring the proper- 
ties of felt or Jello. 

According to LV99, magnetic reconnection of a turbu- 
lent field happens at a rate (see Appendix A) 

Vrec,LV^9~VA{L^.lLY'^Ml for < L, (1) 

where Lx is the horizontal extent of the reconnection 
zone (see Figure [T]), L is the injection scale of turbu- 

^ In special circumstances when the initial configuration of mag- 
netic field contains a magnetic reversal and is laminar, the plasmoid 
reconnection can induce turbulence in 3D (Karamabadi 2012). 
This may be important process for Solar flares (LV99, Lazarian 
& Vishniac 2009, Eyink et al. 2011). However, observations (see 
Larson 1981, Armstrong et al. 1995, Chepurnov & Lazarian 2010, 
Gaensler et al. 2011, see also Elmegreen & Scalo 2006 and Lazar- 
ian 2009 for reviews) testify that interstellar media and molecular 
clouds are turbulent which makes LV99 mechanism directly rele- 
vant. 



lence, Va is the Alfven velocity and Ma = Vl/Va is the 
magnetic Mach number. In Eq. ([1]) it is assumed that 
turbulence is subAlfvenic. For > L, LV99 got the 
expression that differs from Eq. ([1]) by a change of the 
power 1/2 to -1/2. 

The enhancement of the reconneclion rate compared 
to the Sweet-Parker model (see upper plot in Figure [J) 
is achieved by increasing the outflow region A by ac- 
counting for magnetic field wandering. Indeed, in the 
Sweet-Parker model the reconneclion speed is limited by 
the outflow of the matter from the reconneclion layer and 
A is determined by Ohmic diffusion. Therefore, 

Vrec = VA-p- (2) 

is much less than Va since A ^ L^- At the same time A 
becomes comparable with L for V ^ Va provided that 
A is determined by the turbulent field wandering, which 
agrees well with Eq. 

The rate given by Eq. ([T]) is sufficient to allow 
the magnetic field to disentangle during the turnover 
of the eddies. Such eddies are a part of the picture 
of strong Alfvenic turbulence (see Goldreich & Sridhar 
1995, henceforth GS95). The Alfvenic cascade develops 
in compressible MHD turbulence independently from the 
cascade of fast and slow modes, and its properties in com- 
pressible and incompressible turbulence are very similar 
(Cho & Lazarian 2002, 2003, Kowal & Lazarian 2010). 
We use this cascade in its GS95 description for the rest 
of the papei0. For our arguments related to reconnec- 
tion diffusion the compressible motions are of secondary 
importance. Indeed, the LV99 reconneclion model is gov- 
erned by field wandering induced by the Alfvenic cascade 
and turbulent mixing and also depends on the solenoidal 
component of the fiuid. 

The peculiarity of reconneclion diffusion is that it 
requires nearly parallel magnetic field lines to recon- 
nect, while the textbook description of reconneclion fre- 
quently deals with anti-parallel description of magnetic 
field lines. One should understand that the situation 
shown in Figure [1] is just a cross section of the magnetic 
fiuxes depicting the anti-parallel components of magnetic 
field. Generically, in 3D reconnection configurations the 
sheared component of magnetic field is present. 

As we discuss below reconnection diffusion is closely 
connected with the reconnection between adjacent 
Alfvenic eddies (see Figure [3|). The GS95 model predicts 
that Alfvenic turbulence is strongly anisotropic with ed- 
dies elongated along the local direction of magnetic field. 
In other words, the extent of the eddy parallel to mag- 
netic field Z|| is larger than the extent in the perpen- 
dicular direction l±. Mixing motions in MHD turbu- 
lence require that reconnection events in MHD turbu- 
lence should happen through every eddy turnover. This 
is what the LV99 model predicts. Indeed, for small scales 
magnetic field lines are nearly parallel and, when they 
intersect, the pressure gradient is not V2/l\\ but rather 
/l^^)V2, since only the energy of the component of the 
magnetic field that is not shared is available to drive the 
outfiow. On the other hand, the characteristic length 

® Recent studies in Beresnyak (2011) and Beresnyak & Lazarian 
(2010) indicate that the worries that the GS95 description may not 
present a perfect match for numerical simulations are premature. 
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Before reconnection yz view ^ reconnection 




Fig. 2. — Motion of matter in the process of reconnection dif- 
fusion. 3D magnetic flux tubes get into contact and after recon- 
nection plasma streams along magnetic field lines. Left panel: XY 
projection before reconnection, upper panel shows that the flux 
tubes are at angle in X-Z plane. Right panel: after reconnection. 

contraction of a given field line due to reconnection be- 
tween adjacent eddies is l±/l\\- This gives an effective 
ejection rate of Va/1\\- Since the width of the diffusion 
layer over the length Z|| is l±, Eq.([T]) should be replaced 
by Vr ~ Va{1±/1\\)- This provides the reconnection rate 
Va/1\\, which is just the nonlinear cascade rate on the 
scale This ensures self-consistency of critical balance 
for strong Alfvenic turbulence in highly conducting fluids 
(LV99). 

Consider first a toy model illustrating how can plasmas 
move perpendicular to the mean inhomogeneous mag- 
netic field (see Figure [2|). In the figure two magnetic fiux 
tubes with entrained plasmas intersect each other at an 
angle and due to reconnection the identities of magnetic 
field lines change. The process of reconnection requires 
that the magnetic field lines are not parallel. The pro- 
cess is three dimensional with magnetic field component 
reversing its direction in the Z-direction. Thus in XY 
plane we observe nearly parallel magnetic fields, and the 
fact that magnetic field flux tubes cross each other at an 
angle is seen in the subpanel of Figure [H Magnetic re- 
connection feeds upon free energy associated with the re- 
versing magnetic field making magnetic fields smoother. 
For instance, it may be energetically favorable for mag- 
netic flux tubes to reconnect, removing the Z-component 
of the reversed fleld and create a smaller reverse compo- 
nent of the X-field (see the right part of the panel). The 
randomness of the magnetic field in this case decreases, 
which well corresponds to the expectations of the dissi- 
pation process. 

Let us assume that before the reconnection pressure of 
ionized gas Pgas is larger in the left tube (red dots in Fig- 
ure [2) , but the total pressure Pgas + Pmagn is the same 
for the two tubes. In the presence of reconnection, the 
flux tubes get connected and ionized gas streams from 
what used to be the red tube flows to the parts of the 
green tube. In other words, the process of reconnection 
changes the topology of the initial magnetic conflguration 
and connects magnetic fields with different mass loading 
and gas pressures. As a result, ionized gas streams along 
newly formed magnetic field lines to equalize the pres- 
sure along flux tubes. Portions of magnetic flux tubes 
with higher magnetic pressure expand as plasma pressure 
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Fig. 3. — Reconnection diffusion: exchange of flux with entrained 
matter. lUustration of the mixing of matter and magnetic fields due 
to reconnection as two fiux tubes of different eddies interact. Only 
one scale of turbulent motions is shown. In real turbulent cascade 
such interactions proceed at every scale of turbulent motions. 

increases due to the flow of plasma along magnetic field 
lines. The entropy of the system increases with magnetic 
and plasma pressures becoming equal within two newly 
formed flux tubes. This is the diffusion process mediated 
by reconnection that we term "reconnection diffusion" . 

Turbulence is a vital ingredient of reconnection diffu- 
sion. It both induces fast reconnection and transports 
matter and magnetic fields via turbulent motions. Be- 
cause of turbulence, in spite of action of reconnection to 
smooth magnetic field, the magnetic fields cannot relax 
to minimal energy state. Tangled magnetic field always 
emerges and magnetic reconnection inevitably follows. 
The process of entangling magnetic fields by eddies is 
illustrated by Figure [3] where it is shown how two adja- 
cent turbulent eddies can induce magnetic field and mat- 
ter diffusion. There while considering what is happening 
with two flux tubes we see how crossing flux tubes similar 
to those in Figure [2] naturally emerge all the time. The 
eddies of different sizes provide mixing of ionized gas at 
different scales, inducing the diffusion perpendicular to 
the mean magnetic field that is not limited by the speed 
of ionized gas moving along magnetic field lines. 

Indeed, if the densities of plasma along magnetic field 
lines are different in the two flux tubes, the reconnec- 
tion in Figure [3] creates new flux tubes with columns 
of entrained dense and rarefied ionized gas. The situ- 
ation is similar to that in Figure [D and the gas moves 
along magnetic fields equalizing the pressure within the 
newly formed flux tubes. As a result, eddies with initially 
different plasma pressure exchange matter and equalize 
plasma pressure. This process presents the diffusion of 
ionized gas perpendicular to the mean magnetic field and 
the rate of this diffusion does not depend on the gas ion- 
ization. The process happens at all scales of self-similar 
eddies constituting the turbulent cascade. The alterna- 
tive microscopic picture of the diffusion of magnetic field 
lines is provided in Appendix B. 

In reality, for turbulence with the extended inertial 
range, the shredding of the columns of ionized gas with 
different density proceeds at all turbulence scales mak- 
ing speed of plasma motion along magnetic field lines 
irrelevant for the diffusion. For the case of strong tur- 
bulence the diffusion of matter and magnetic field over 



large scales > L is given by Eq. ^ . 

Mixing motions of magnetic eddies induce turbulent 
diffusivity with the diffusion coefficient (Lazarian 2006): 

^strong ^ 

LVlMI (3) 

where the reduction compared to the usual hydrody- 
namic turbulent diffusivity coefficient, i.e. 

Khydro ^ LVl , (4) 

arises from a relative inefficiency of mixing induced by 
weak turbulence. The latter is induced on the injec- 
tion scale for isotropic subAlfvenic turbulent driving, i.e. 
Ma < 1. For superAlfvenic turbulence the estimate 
given by Eq. Q is appropriate. 

Naturally, the estimate of diffusivity given by Eq. ([3]) 
is applicable if the scales that we deal with are larger 
than the injection scale of strong MHD turbulence. The 
latter is also the scale for the transition from weak to 
strong turbulence (LV99, Lazarian 2006): 

krans ^ L{Vl/Va)^ = LM^. (5) 

If the scales that we deal with are smaller than the scale 
given by Eq. ([SJ, then we deal with the accelerating 
process of turbulent diffusion. Indeed, the separation 
between two particles dl{t)/dt ^ v{l) for Kolmogorov 
turbulence is atl^^^ , where a is proportional to a cube- 
root of the energy cascading rate, i.e. at ~ Vl/L for 
turbulence injected with superAlvenic velocity Vl at the 
scale L. The solution of this equation is 

l{t) = [ll^' + at{t-to)]'^', (6) 

which provides Richardson diffusion or ~ t^. The ac- 
celerating character of the process is easy to understand 
physically. Indeed, the larger the separation between the 
particles, the faster the eddies that carry the particles 
apart. The character of the Richardson diffusion in terms 
of field line separatiorQ is illustrated by the left panel of 
Figure [TT] in Appendix C. 

When we deal with the turbulence within a cloud, it 
is plausible to associate the injection scale of the turbu- 
lence with the size of the cloud. Indeed, although the 
clouds may be a part of the compressible turbulence cas- 
cade originating at large scales, the disparity between 
the Alfven velocities in dense and rarefied gas mitigates 
the cross-talk between the eddies inside and outside the 
cloud. 

3. OBSERVATIONS OF MAGNETIC FIELDS IN 
MOLECULAR CLOUDS 

Star formation is known to be a rather inefficient pro- 
cess. Indeed, the interstellar mass of the Milky Way is 
Mmw ~ 10^ Mq. For the typical density of the gas of 50 
cm-3 the free fall time is t// « (37r/32Gp)i/2 ~ 5 x 10^ 
years, which provides a "natural" star formation rate 
MMw/Tff of 200 Mq per year. At the same time the 
measured star formation rate is « 1.3 Mq per year (Mur- 
ray & Rahman 2010). It is because of this inefficiency 
that we still have large amounts of interstellar media in 

^ The description of the Richardson diffusion of magnetic field 
lines was provided in LV99. 
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spiral galaxies. A traditional way of explaining this incfB- 
ciency is to appeal to magnetic forces preventing gravita- 
tional collapse (see Mouschovias & Spitzer 1976). Indeed, 
it is possible to show that if the ratio of the magnetic flux 
to mass is larger than the critical one, 

($/M)erit « 1.8 X 10"^ gauss cm^ g-^, (7) 

the magnetic field prevents cloud collapse (see Draine 
2011) and star formation should come to a halt. If am- 
bipolar diffusion mediates star formation, the time scales 
are approximately right, which was a strong argument in 
favor of the ambipolar diffusion paradigm. 

A crucial parameter for star formation theory is tlicac;- 
fore the ratio of mass to magnetic flux, M/$, which pa- 
rameterizes the relative importance of gravity and mag- 
netic pressure. Strictly speaking, one should include all 
of the mass in a magnetic flux tube to calculate this ratio. 
However, what is important for star formation is the local 
ratio for a self-gravitating cloud. Other mass located far 
away in the flux tube will not affect the balance between 
gravity and magnetic pressure for the cloud. Observa- 
tionally, one measures this local M/<^ by detecting the 
Zecman effect in spectral lines of a species (such as OH 
or CN) that sample the densities and spatial regions of 
the cloud. Within a telescope beam of radius r the mag- 
netic flux being sampled is $ oc wr'^B, while the mass 
of the cloud within that beam is M oc wr'^N, where A'' 
is the column density. Hence, the local M/$ oc N/B. 
Since both B and A'' can be inferred from the Stokes 
V and / spectra, one infers the local that is, the 

M/$ in the region of space and velocity sampled by the 
observed species. This is the appropriate parameter for 
considering the balance between gravity and magnetic 
pressure in a cloud. This local M/$ can change with- 
out magnetic diffusion if there are flows of matter onto 
(or from) the cloud from distant regions of the flux tube. 
Such flows provide a mechanism for building up a local 
self-gravitating mass (a cloud) without changing the lo- 
cal magnetic field strength. M/$ could then vary from 
one place to another within a cloud. One issue is whether 
such flows can operate over the requisite distances with- 
out turbulence and instabilities disrupting the flow. A 
general discussion of such issues are beyond the scope 
of the present paper. However in Appendix D we pro- 
vide analytical arguments showing that the accumulation 
along the magnetic field lines is unlikely in observations 
by Crutcher et al. (2011). Here we consider an alterna- 
tive diffusion process to ambipolar diffusion for changing 
local M/$. 

The ambipolar diffusion theory of star formation is well 
developed, with several specific predictions. Crutcher 
(2012) reviewed the current state of observations of mag- 
netic fields in molecular clouds and compared the obser- 
vations with the predictions. Here we briefiy summarize 
four observational tests of the theory. 

First, ambipolar diffusion predicts and requires that 
M/$ increase from the envelope to core region of a 
cloud. Crutcher, Hakobian & Troland (2009) directly 
tested this prediction. Their experiment measured TZ = 
{M/^)core/{M/^)enyeiope- The idealized ambipolar dif- 
fusion theory of core formation requires TZ to be approx- 
imately equal to the inverse of the original (before evo- 
lution driven by ambipolar diffusion) subcritical M/$, 
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Fig. 4. — HI, OH and CN Zeeman measurements of line of sight 
component of magnetic field versus the total column density Nh 
of atomic and molecular hydrogen. The straight line is the critical 
mass/flux M/^ = 3.8 X Njj/B with subcritical clouds above the 
line. The two molecular clouds near N 10^^ have OH maser con- 
tamination in their spectra and the corresponding Zeeman results 
arc not reliable. From Crutcher 2012. 

or TZ > 1. The test used OH Zeeman observations to 
determine the line-of-sight magnetic field Bios and the 
column density Nqh in each of four dark cloud cores 
and envelopes. The central assumption of this test was 
that the angle 6 between the line of sight and the mag- 
netic field was approximately the same in core and en- 
velope. Then the measured Bios = Btot cos 9 would have 
the same 6 in the numerator and denominator of TZ, 
and the dependence on the unknown 9 would disappear. 

Therefore, TZ = {Nqh I Bios)core/ {Nqh / Bios)envelope, all 

directly measurable quantities. 

Because cloud cores with the strongest measured Bios 
were selected for the test, it is likely that Bios ~ Btot 
with ^ ~ 0, so cos^ ss 1 and small changes in 9 between 
core and envelope would not affect TZ very much. Sup- 
porting this assumption of the experiment is the fact that 
published models of core formation driven by ambipolar 
diffusion have strong, regular magnetic field morphology 
such that the unknown angle 9 between B and the line 
of sight is approximately the same in core and envelope 
regions. Observationally, field morphologies in cores and 
envelopes are found to be correlated in direction (Li et 
al. 2009). 

With measurements of Bios and N{OH) toward each 
core and envelope, it was possible to calculate TZ and its 
uncertainty toward each cloud. Because Bioa was not 
detected (at the 3ct level) toward any of the four cloud 
envelopes, in no case was the calculated value of TZ signif- 
icantly (3(t) different from 0. But the astrophysical qties- 
tion is not whether TZ is significantly different from 0, it 
is whether the TZ were significantly different from 1 , since 
the ambipolar diffusion thcorv requires that TZ > 1. The 
results for the four clouds were: TZ{LIU^C0) = 0.02 ± 
0.36,7^(B217-2) = 0.15 ± 0.43, 7^(L1544) = 0.42±0.46, 
and TZ{B1) = 0.41 ±0.20. Hence, in all four cases, the ex- 
periment found TZ <1, not 7?. > 1 as the theory predicts. 
In two of the four cases, 7^ < 1 at near a Scr level. 

Proponents of the ambipolar diffusion theory of star 
formation strongly attacked the Crutcher, Hakobian & 
Troland (2009) results and conclusions (Mouschovias & 
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Tassis 2009, Mouschovias & Tassis 2010). The four 
strongest objections and the counters to those objectfons 
are addressed below. 

(1) Envelope regions around the core positions are not 
uniform, so the four envelope results should not be com- 
bined to improve the signal-to-noise ratio. For each enve- 
lope position considered separately, the sensitivity to the 
differential M/^ is reduced and the disagreement with 
ambipolar diffusion less significant. However, at all four 
envelope positions in all four clouds, the column density 
is smaller than that at each central position, so envelope 
positions are being measured, just not the identical enve- 
lope column density at each position. The mean envelope 
column density is about 1/2 that of each core position. 
Hence, by obtained a mean Af/^ for each envelope, one 
is measuring mean envelope properties that are entirely 
suitable for comparison with the central core properties 
for each cloud. This procedure is similar to the stacking 
operation used in galaxy studies to infer mean properties 
around the center of galaxies, and is an entirely accepted 
and appropriate technique for obtained mean properties. 

(2) Since Bios was not detected in the envelopes only 
upper limits should be considered. However, using the 
obtained value and its uncertainty to compare with the 
prediction of theory is equally valid. If theory predicts 1 
and observations give 0.5 ± 0.25, it is valid to state that 
the experiment differs from the prediction by 2a, which 
is the approach we take. Their preference for stating that 
the data agree with the theory at the 2a level is not a 
different statement. In no case do Crutcher, Hakobian 
& Troland (2009) claim a result > 3cr, but merely state 
that the data for all four clouds do not correspond to the 
result predicted by theory. Moreover, the assumption of 
the theory that clouds initially are subcritical implies not 
that TZ ~ 1 but that 7?. > 1, so the significance of the 
individual observational results for ambipolar diffusion is 
larger than given by comparison with TZ = 1. 

(3) The motion of cores through surrounding more dif- 
fuse gas could lead to B in cores and their envelopes not 
being essentially parallel. However, while this is certainly 
possible, it would require that all four cores move nearly 
in the plane of the sky so the envelope fields would be 
dragged into the plane of the sky. Also, as mentioned 
above, evidence from linear polarization mapping that 
directions in cores and surrounding gas are strongly cor- 
related (Li et al. 2009) argues against this morphology. 

(4) Multiplication of the probabilities of the data for 
each cloud being consistent with ambipolar diffusion to 
assess agreement with theory is invalid because the four 
clouds are measurements of different cases, so probabil- 
ities cannot be combined. However, while the observa- 
tional results for each cloud individually can to argued to 
be (marginally) consistent with the prediction of ambipo- 
lar diffusion, one must also consider that a sample of four 
clouds was observed. If in fact 7?. > 1 in all four clouds, 
one would expect observational noise to sometimes pro- 
duce an observed TZ greater than the "true" value, not 
always smaller. While the sample size of four clouds is 
not large, the observations reported by Crutcher, Hako- 
bian & Troland (2009) certainly suggest a problem with 
the idea that ambipolar diffusion was responsible for the 
formation of cores in all four clouds. If the probability 
for one cloud is pi, it is correct to calculate piP2P3Pi 
as the probability that all four were formed by ambipo- 



lar diffusion. Hence, the probability stated by Crutcher, 
Hakobian & Troland (2009) that all four of the clouds 
have 72. > 1 is 3 X 10"'' is statistically correct. 

Second, although the M/<I> ratio is an input parameter 
to the ambipolar diffusion theory of star formation, for 
the theory to be relevant clouds must initially be sub- 
critical, with ambipolar diffusion leading to an increase 
of mass in cores until cores become supercritical and col- 
lapse. The physics of the theory would be perfectly valid 
if subcritical clouds did not generally exist in nature, but 
then ambipolar diffusion would not be the physics driving 
star formation. Zeeman data make it possible to assess 
whether relatively low-density self-gravitating molecular 
clouds are generally subcritical. Figure |4] from Crutcher 
(2012) shows Zeeman observations of Bios versus Nh- 
Although an individual Zeeman measurement gives Bios 
(hence only an upper limit to A//$), the upper envelope 
of the Bios in Figure [4] defines the total field strength 
Btot at each Nh, since for some fraction of the clouds B 
should point approximately along the line of sight. For 
Nh ^ 10^^ cm~^, the M/<I> are subcritical; these clouds 
are lower density H I clouds. For Nh ^ 10^^ cm~^, 
the M/$ are overwhelmingly supercritical; these clouds 
are higher density molecular clouds and cores. Hence, 
the data appear consistent with the theory, with neu- 
trals gravitationally contracting while leaving the mag- 
netic flux behind and hence increasing M/$ in the higher 
density molecular gas. However, the H I clouds are 
in approximate pressure equilibrium with the warm in- 
terstellar medium and are not self-gravitating, so they 
could not gravitationally collapse through the magnetic 
field to form supercritical cores. For Nh 10^^ cm~^ 
most of the points in Figure [4] are molecular clouds that 
are generally self-gravitating, so this should be the re- 
gion of transition from subcritical to supercritical clouds. 
Yet there are zero definite cases of subcritical clouds for 
Nh > 10^^ cm~^! The two points which seem to be 
significantly above the critical line are both from OH 
in absorption toward H II regions, with polarized maser 
emission within the absorption lines producing inconsis- 
tencies and uncertainties in the fit for Bios that were not 
reflected in the formal uncertainties. Hence, there are 
no clear observational cases of self-gravitating, subcriti- 
cal clouds, which must exist in the ambipolar diffusion 
theory of star formation. 

The third and fourth tests come from a statistical 
analysis of Zeeman data. For ambipolar diffusion to 
be the principal driver of star formation, magnetic field 
strengths in molecular clouds should be strong, so that 
initially clouds are subcritical. Ambipolar diffusion does 
not decrease field strengths, it has neutrals collapsing 
through magnetic fields essentially leaving the field in 
place until the core becomes supercritical. The field 
strength in the core would then increase with density. 
For & B IX scaling parameterization, k would increase 
from near to 0.5 as the collapse proceeded. Hence, 
one would expect the total magnetic field strengths in 
molecular clouds to be high so M /$ would be near criti- 
cal, and the field to be sufficiently strong that it controls 
the power law exponent k to be within the above range. 
Crutcher et al. (2010) carried out a Bayesian statistical 
analysis of Zeeman data in order to infer the PDF of Btot 
from the observed Bios and compared results with these 
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Fig. 5. — H I, OH, and CN Zeeman measurements of the magni- 
tude of Bjos versus n/^ . The dashed hne segments show the most 
probable model from the Bayesian analysis for the maximum total 
magnetic field strength at each density. The total field strength of 
an individual cloud lies between this maximum and essentially 0, 
with a fiat distribution between and the maximum; from Crutcher 
et al. (2010b). 

final two predictions of tfie ambipolar diffusion tlieory. 
Figure [S] sliows tlie Zeeman measurements of Bios versus 
hh used in tfie Bayesian analysis. 

The third prediction of ambipolar diffusion is that most 
molecular clouds should be subcritical, except for those 
that have achieved supercritical cores and are collaps- 
ing; no very weak magnetic fields should be found. The 
most probable PDF of Btot was found to be a flat func- 
tion; that is, for a set of clouds at any density uh the 
Btot were uniformly distributed between very small val- 
ues and the maximum Btot found for that nn (the dashed 
line shown in Figure [S] The mean M/ $ for molecular 
clouds was found to be supercritical by ~ 3. Hence, as 
noted above, no subcritical molecular clouds have been 
shown to exist, so all > 1 (supercritical), and many 

molecular clouds have very small Btot and hence are very 
supercritical. 

Fourth, the Bayesian analysis found that k w 0.65 ± 
0.05 for the B oc scaling of total field strength with 
density. The ambipolar diffusion theory predicts k < 0.5. 
A fc = 2/3 value would be produced by contracting clouds 
with weak magnetic fields, such that magnetic pressure 
would not significantly affect the geometry of the collapse 
(Mestel 1966). Such weak fields are consistent with the 
result for the PDF of total field strengths just discussed. 

There are thus four discrepancies between the predic- 
tions of ambipolar diffusion theory and observations of 
magnetic fields in regions of star formation for the exist- 
ing paradigm of star formation - near flux freezing with 
ambipolar diffusion driving contraction of cloud cores. 
In the next two sections we will consider whether re- 
connection diffusion theory is consistent with these four 
observational results. 

4. SIMPLE MODEL FOR RECONNECTION DIFFUSION 

We first consider the first test of the above section. 
We will directly compare a reconnection diffusion model 
with the results of the Crutcher, Hakobian & Troland 
(2009) observations that M/$ decreases rather than in- 
creases between envelopes and cores of the four clouds 
they studied. The relevant mean parameters for the four 
cores are Btotai ~ 30 /zG, nn ~ lO"* cm~'^, and radius 



w 0.2 pc, and for the envelopes Btotai ^5—10 /iG, 
riH ~ 10'^ cm~'^, and radius ~ 1 pc. If these clouds 
formed by contraction of a typical diffuse H I cloud, the 
parameters would be Btotai ~ 5 /iG, uh ~ 60 cm~^, and 
radius ~ 2.5 pc. 

As we mentioned earlier, in the presence of grav- 
ity reconnection diffusion allows density to concentrate 
towards the gravitational center, while magnetic field 
leaves the gravitational potential. The reconnection dif- 
fusion rate and thus the rate of magnetic field removal 
is expected to be proportional to the level of turbulence. 
Simulations that illustrate this effect are presented in 
Figure [51 It is known that the velocities in cloud en- 
velopes are larger than the velocities in the cloud cores 
(see Taffalla et al. 1998). In addition, the scale of tur- 
bulence involved in reconnection diffusion is also larger. 
Thus we expect the diffusion coefficient for the envelope 
(see Eq. |3]) to be larger compared to the core coefficient, 
and it is natural that faster removal of magnetic field 
happens from the envelope than from the core. In fact, 
we also expect the uniformly compressed magnetic field 
to be smaller at the edges of the cloud even for uniform 
diffusivity. These theoretical expectations agree with the 
results in Crutcher, Hakobian & Troland (2009). The ob- 
servational results disagree with the ambipolar diffusion 
paradigm as the rate of ambipolar diffusion decreases to- 
wards edges of clouds, and within clouds the ambipolar 
diffusion of magnetic fields is negligible for media with 
high degree of ionization. On the contrary, reconnection 
diffusion does not depend on the degree of ionization, 
and it can successfully proceed in highly ionized gas. 

In what follows, we illustrate the effects of reconnec- 
tion diffusion using simplified calculations of magnetic 
diffusivity for conditions motivated by the observations 
in Crutcher et al. (2009). 

4.1. Solutions of the model: the radial profile produced 
by magnetic reconnection 

Magnetic fields are weightless. Therefore we can get 
an insight in the physics of the reconnection diffusion by 
modeling the diffusion of magnetic field from the place 
of its original concentration. For illustrative purposes 
we shall limit our discussion to the one dimensional dif- 
fusion: 

dB d^B „ 

= 0. (8) 



dt 



This equation can be solved analytically for a constant 
diffusion coefficient k and an initial distribution of mag- 
netic field 

B{x,t = 0) = (l){x). (9) 
The solution of the equation is given by 

cj^ix') 



B{x,t) 



-oo 2-\/ IT K,t 



■i^t dx , 



(10) 



which can be easily solved for the case of a uniformly 
compressed fluid 



^ ( Bo, for|a:|<i? 
\ BisM, for > R, 



(11) 



where Bq denotes the initial value of the magnetic field 
after the compression, and Bism the magnitude of the 
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Fig. 6. — Left panel. Schematic of the cloud core and envelope in Crutcher et al. (2009). Both scale of the turbulence and velocity (and 
thus reconnection diffusion rate) are larger in the envelope than in the core. Right panel. Results of 3D MHD simulations (modified from 
Santos-Lima et al. 2010) The removal of magnetic field via reconnection diffusion increases with the increase of the turbulence level. 



field outside. The solution can be obtained considering 
a constant boundary condition (i.e. B(x,t) — Bism, for 
\x\ > R), and it is given by the error function 

B{x, t) = SisM + (So - Sism) X 
' R-x" 



Erf 



Erf 



R 



(12) 



Let us estimate the reconnection diffusion coefficients 
based on the observations of Crutcher et al. (2009). We 
have considered a cloud core with size (radius) £ = 0.2 pc, 
within an envelope of radius L — \.Q pc. For a magnetic 
field of B = 10 /iG in the envelope, and a hydrogen 
number density uh — 10"^ cm^"^, the Alfven speed is 



^^,cnv — 



B, 



,cnv 



0.6 km s 



(13) 



for fi — 1.3 (neutral gas with solar abundances). For the 
core we assumed i? = 35 ^G, and nn = 10"* cm"'^, which 
results in a similar VA.coro ~ 0.67 km s^^. 

For the velocity we use the mean dispersion obtained 
from the OH lines in the clouds observed by Crutcher 
et al. (2009). They measure radial velocities with a 
mean full width at half maximum of 0.99 km and 0.76 
km s~^, for the envelope and core, respectively. These 
correspond (dividing by 2V2 In 2 « 2.35) to a velocity 
dispersion of 0.42 km s^^ and 0.32 km s^^, respectively. 
Assuming that the motions are isotropic, the magnitude 
of the total velocity dispersion can be obtained by mul- 
tiplying those values by ^/3, yielding Vl ~ 0.73 km s~^ 
in the envelope, and Vg = 0.55 km s~^ in the core. 

The resulting Alfvenic Mach numbers are Ma.cuv ^ 1-2 
for the envelope, which falls into the strong turbulence 
regime, and Ma^covc ~ 0.8, corresponding to the weak 
turbulence case. 

Let us illustrate the diffusion of magnetic field due to 
reconnection by taking two discrete values for the n co- 
efficient, one for the core Kcorc (for weak turbulence) and 
one for its envelope Kcnv (for strong turbulence), such 
that: 

Kcore =V^^^Mi ,„,,w 1.86 X lO^^ cm^ s-\ (14) 



= VlL 



2.25 X 10 



23 cm^ s-i. 



(15) 



Which gives a ratio of i^cmj^corc ^ 12; thus the rate 
of removal of magnetic field through turbulent reconnec- 
tion is roughly twelve times faster with Hcnv, with a dif- 

T 2 

fusion timescale of tdiff,cnv = -j: — ~ 1.3 Myr. In Figure 
|7| we illustrate the solution of Eq. (IT^ for the initial 
conditions given in Eq. (fTTj) . with R — L = 1.0 pc, 
Bo = Bcorc = 35 /xG, and -Bism = 6 fiG. In the figure 
we see how the magnetic field diffuses from the envelope 
towards the core, producing already a profile that is con- 
sistent with the observations (i.e. a larger field in the 
core than in the envelope). Of course, this is an overly 
simplified model, as the diffusion coefficient is not con- 
stant throughout the region. It gives, however, an insight 
of the timescales involved in either case. 

We also search numerically for a solution assuming that 
the diffusion coefficient k, scales as P, a > 0, which 
represents the increase of the diffusion efficiency with 
the scales, which is expected from turbulent diffusivities. 
We assume for simplicity that a = 4/3 corresponding to 
the GS95 picture of turbulence. The scale dependence is 
K oc ^^Z"^, where the proportionality constant was chosen 
so that K = Kcorc for x — 0, and k = Hcnv for x = L. 

The solutions (numerical) with a scale dependent co- 
efficient are presented in Figure |8| for two different initial 
conditions. In the top panel, the same field used in Fig- 
ure |7| is used. In the bottom panel we have used: 



35 nG, 
6 AiG, 



for 
for 



< ^ 
> i. 



(16) 



In this solutions we can see the effect of a more efficient 
removal of magnetic field in the envelope than in the core. 
In the case in which the initial field was concentrated 
at the core (Figure Hb) it is also evident how some of 
the field from the core diffuses to the envelope. In that 
case and after ~ 200 kyr, the magnitude, both for the 
core and the envelope are comparable to the observations 
of Crutcher et al. (2009). An important point is that, 
regardless of the initial condition, the profiles preserve 
an increasing magnetic field towards the center, with a 
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Fig. 7. — DifTussion of t he m agnetic field with a constant diffu- 
sion coefficient (Equation If 21 1. fn the top panel (a) the solution 
is obtained with the turbulent diffussion coefficient Kcnv, and in 
the bottom panel (b) the same solution, but with the coefficient 
in the core (kcofg) is presented. The different lines, as indicated 
by the color coding in the legend, correspond to different times at 
intervals of 65 kyr, from t = to t = 455 kyr. 

larger strength ratio between center and envelope than 
ambipolar diffusion would produce. 

5. OTHER OBSERVATIONAL TESTS 

5.1. No subcritical molecular clouds 

The second test discussed in §3 is the fact that there 
is no observational evidence for the subcritical self- 
gravitating molecular clouds predicted by the ambipo- 
lar diffusion theory to be observable among the sam- 
ple of several dozen dark clouds (densities 10'^~'* cm^'^) 
for which OH Zeeman observations are available. Such 
clouds are the initial condition for the ambipolar diffu- 
sion theory to be relevant to star formation. On the 
other hand, this observation is a natural consequence of 
reconnection diffusion, which acts in low-density clouds 
where turbulence is strong to prevent an increase of mag- 
netic field strength with density. In low-density clouds 
the kinetic and magnetic energies are approximately in 
equipartition and larger than the thermal energy (Heiles 
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Fig. 8. — Diffusion of the magnetic field with a scale dependent 
diffusion coefficient, k oc ^^Z^. The initial condition for the mag- 
netic field is represented by the black line, a single step in panel 
(a), and a two step in panel (b), see text for details. The different 
lines, as indicated by the color coding in the legend, correspond to 
different times at intervals of 32.5 kyr, from t = to i = 227.5 kyr. 



& Troland 2005). From the point of view of our anal- 
ysis, this is a signature of the transAlfvenic supersonic 
turbulence indicating that the turbulence should be effi- 
cient in driving reconnection diffusion. We note, how- 
ever, that an alternative picture for building up self- 
gravitating clouds without significantly increasing mag- 
netic field strengths is the collection of matter along mag- 
netic field lines. In a simplistic picture of this process, 
density would increase while the magnetic field strength 
would be unchanged. Such a process would produce thin 
cylinder or disk-like clouds with the magnetic field along 
the minor axis. However, there is not universal observa- 
tion support for this morphology. Goodman et al. (1990) 
and Chapman et al. (2011) found clouds with minor axes 
parallel, perpendicular, and intermediate angles with re- 
spect to the surrounding magnetic field. A statistical 
analysis by Tassis et al. (2009) found that the mag- 
netic fields were oriented close to 24°) the minor axis 
but could not reject alternative configurations. No doubt 
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matter can accumulate more along than perpendicular to 
magnetic fields. But the scales for such one-dimensional 
collection to form GMCs are of hundreds of parsecs (see 
Vazquez-Semadeni et al. 2011, Appendix D.). Field wan- 
dering induced by turbulence (see Appendix D) would 
tend to interfere with a smooth flow over such distances. 
We believe that reconnection diffusion is an intrinsic part 
of the simulations of the turbulent interstellar medium, 
and it should be accounted for in interpreting the results 
of turbulent numerical simulations. 

5.2. Clouds with very weak fields 

The third observational result that ambipolar diffu- 
sion theory has difficulty with is the fiat PDF of total 
field strengths (Crutcher et al. 2010); that is, at a given 
density field strengths range from very small values cor- 
responding to a highly supercritical M/$ up to values 
that correspond to an approximately critical M/$. The 
molecular clouds with very small field strengths could 
arise from the action of reconnection diffusion over time, 
with older clouds having field strengths approaching that 
of the surrounding H I medium. While ambipolar diffu- 
sion theory does not prohibit the observed highly super- 
critical molecular clouds at relatively low densities (10'^^^ 
cm~^), for the theory to be relevant most clouds must 
start as subcritical. 

5.3. Scaling of field strength with density 

The fourth observational result, that field strengths 
scale as density to the 2/3rds power at high densities 
(Figure [SJ, was predicted by Mestel (1966) for contract- 
ing clouds with magnetic fields too weak to significantly 
affect the collapse. We believe that the necessary low 
field strengths are due to reconnection diffusion being ef- 
ficient at the lower densities before collapse. Once rapid 
collapse begins, reconnection diffusion removes magnetic 
field from the contracting clouds more slowly than con- 
traction increases it. If the reconnection diffusion time 

trec.diff = I'^/n , (17) 

where / is the cloud size, is less than the free fall time 



i.e. t 



tff = V37r/(32Gp) , 



(18) 



reconnection diffusion should be efficient in keeping the 
magnetic field inside the cloud close to the value of the 
magnetic field outside the cloud (see discussion in §5.2). 

If this is true, for the relevant scales reconnection diffu- 
sion must be faster than free fall. Consider the scales at 
which trec,di f f (scc Eq. (ITTI) ) is less than the free fall time 
tff (see Eq. (jlSp ) and compare our expectations with 
observations. If for the case of subAlfvenic turbulence, 
one uses the scale I of the cloud and the corresponding 
velocity 

vi=V,njil/L,r.jy/^ (19) 

where Li„j and Vinj are the injection scale and velocity, 
respectively, then the Kolmogorov scaling of velocities 
given by Eq. (|T9|) corresponds to the predictions of the 
GS95 model for scales less than krans = Linj{VL/VA)^ 
(see Lazarian 2006). Thus one can get the condition 
that reconnection diffusion is faster than cloud free fall. 



rec.dif f 



llower ^ 



< t 



ff 



if 



Linj (Gp) 



1/2 



Va 



(20) 

In the case that magnetic field stays roughly constant as 
in Figure [U the Alfven velocity Va in Eq. ([20| changes 
with density as p"^!"^ , and the requirement on the scale 
of the cloud can be rewritten 



II 



> 



1 



7^663/2 y-12 



(21) 



which provides the lower constraint on the size of the 
cloud of given density for which reconnection diffusion is 
efficient. 

If turbulence is superAlfvenic at the injection scale 
Linj, at small scales it gets subAlfvenic. If one approxi- 
mates superAlfvenic motions as well by Kolmogorov type 
scaling (given by Eq. (fT9)) ). it is possible to find the con- 
straint on I for the motions to be subAlfvenic at the scale: 



''Sub 



(22) 



For scales larger than this, the turbulence is super- 
Alfvenic and the reconnection diffusion is described by 
the turbulent diffusivity coefficient ^ Ivi, and the cri- 
terion for the reconnection diffusion to dominate at the 
scale I is 



lu 



pper 



< 



37t\ 
32 / 



3/4 x.3/2 

inj 



1 



r 1/2 

''inj 



[Gpf/^ 



(23) 



Within the range of [kowenlupper], the reconnection dif- 
fusion induced by turbulence removes magnetic flux from 
clouds before they collapse. Beyond this range the gravi- 
tational compression of matter will increase the strength 
of the magnetic field in agreement with the tendency 
shown in the observational data in Figure [S] 

In Figure [H] we illustrate the region delimited by 
[hower, lupper] for a fcw selcctcd models of turbulence in- 
jection. In particular, we consider Model 1 with L^j — 
30 pc, with Vl = 10 km s~\ Model 2 with L.„y = 100 pc, 
with Vl = 10 km s^^ and Model 3 with L^j = 30 pc, 
with Vl = 20 km s^^. In spite of the decades of studies, 
the characteristics of turbulence are not well determined. 
We believe that Model 1 may be closer to reality in terms 
of representation of the parameters of turbulence in the 
galactic disk. It conservatively adopts the scale which is 
smaller than the one obtained in Chepurnov et al. (2010) 
for atomic HI at high latitudes. In other words, it takes 
into account that the scale of turbulent eddies may be 
smaller in the disk compared to high latitudes. In two 
panels of Figure[n]the results are illustrated for two values 
of magnetization, namely, for B — 6 pG and B = 30 pG. 
Naturally, for Model 1 the calculations are valid only up 
to the scales less than injection scale, which is shown by 
horizontal lines in the upper parts of the panelf0. Figure 
[9] also shows the regions of subAlfvenic and superAlfvenic 
turbulence. The boundary between the latter regions is 
given by Eq. 

* Therefore for Model 1 the hnes do not extend beyond the hor- 
izontal Une. 
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Fig. 9. — Illustration of the parameter space for which reconnection diffusion is faster than the gravitational collapse. The plots corresponds 
to three different models of turbulent injection described in the text. Lines corresponding to iuppei-, hower and the Jeans scale Rj are 
shown. The horizontal lines corresponds to the injection scale of 30 pc of Model 1. Left panel: Case of magnetic field o{ B = 6fj,G. Right 
panel: Case of iJ = 30^G magneti field. 



For illustrative purposes we also show on the plots the 
Jeans scale 

where iiniH, is the molecular mass of the gas, which is 
2 in our case, and T is the gas temperature T, which is 
20K for our case. The corresponding line in Figure [9] in 
most cases is above the lines of hower and therefore no 
free fall happens for such small lengths. That is why we 
focus in our discussion on lupper boundary. 

According to Figure [5] the transition to the regime in 
which reconnection diffusion is slower than free fall hap- 
pens at roughly n w 10* cni~^. The critical cloud scale 
corresponding to the dominance of the free fall is given 
by lupper boundary. For Model 1 in Figure [9] the scale 
is around several parsecs. According to Eq. (1^ this 

scales depends as L~^-'^ , which results in small varia- 
tions within the range of possible injection scales. We 
see that within the range of parameter variations illus- 
trated by Figure |9] the critical cloud scale does not vary 
significantly. Thus for our models the column density 
of the relevant clouds is expected to be of the order of 
Nh ~ 10^^ cm~^. While this should be treated only as 
a very rough estimate due to the nature of our simpli- 
fied treatment, it is encouraging that this column density 
corresponds to the observational data in Figure ID 

In the scenario that we advocate here the clouds that 
undergo further collapse are supercritical and therefore 
magnetic fields cannot stop the free fall. The intensity of 
the magnetic field is increasing, however. In reality, the 
approximation of the cloud collapse with i is an ideal- 
ization. Clouds being steered by turbulence get virialized 
and free-falling clouds do not exist. Therefore, in agree- 
ment with Figure 21 we expect to see dense clouds which 
existed long enough for reconnection diffusion to remove 
the excess of compressed magnetic flux. As the density 
increases, gravity becomes more important and the ten- 
dency of increasing magnetic field with density gets more 
prominent in Figured 



Our analysis above is of illustrative nature and we ac- 
cept the existing uncertainties related to the injection 
scale and velocity as well as influence of the compress- 
ible motions and the effect of compressions on turbulence 
(see Robertson & Goldreich 2012). Nevertheless, our re- 
sults indicate the existence of the domain [hower, lupper] 
within which the reconnection diffusion should be faster 
than the free fall. Outside the domain, cloud collapse can 
happen more rapidly than reconnection diffusion, drag- 
ging the magnetic fleld with it. This may be responsible 
for the 2 /3rds power-law dependence of fleld strength on 
density that we observe (Figure [5]) for densities larger 
than n ^ 10* cm~^. 

6. FURTHER IMPLICATIONS OF RECONNECTION 
DIFFUSION 

6.1. Larson laws 

We note that the observational results reviewed by 
Crutcher (2012) on the approximate independence of 
magnetic fleld strength on the cloud column density up 
to Nh ~ 10^^ cm~^ (and the reconnection diffusion con- 
cept that accounts for these results) can be related to 
the empirical Larson relations (Larson 1981) obtained 
for interstellar turbulence. Larson (1981) found that the 
velocity dispersion is proportional to the square root of 
the cloud size, i.e. 

, (25) 

and that the 3D density of a cloud is inversely propor- 
tional to cloud size, i.e. 

p^R^^ . (26) 

For instance, one can assume a rough equality between 
the kinetic energy and magnetic energy 

£ - ^^^^ ' 

which is true for transAlfvenic turbulence where motions 
are associated with non-linear Alfven modes (Goldreich 
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& Sridhar 19950 It is also natural to accept the virial- 
ization of a cloud 



The Eqs. (Ull) (HH), if combined with the simplest esti- 
mate of cloud mass M ~ pB? , yield 

a. ^ B^'^R^'^ (29) 

and 

pr^BR-\ (30) 

These reproduce the Larson (1981) relations given by 
Eqs. (|25|) and (|26p if, in agreement with observation^^ 
(see Figures [5] and H} , the reconnection diffusion keeps 
magnetic field uncorrelated with density. 

For cloud cores where reconnection diffusion is not fast 
enough to remove magnetic field on the time scale of the 
dynamic collapse (see Tafalla et al. 1998, Reiter et al. 
2011), the Larson relations fail, in agreement with ob- 
servations and simulations (see Nakamura & Li 2011). 
Quantitatively, the criterion for the reconnection diffu- 
sion to be able to remove the magnetic field from the 
collapsing cloud of size Idoud is Vinfaii < n/ldoud, where 
K is one of the diffusion coefficients (see Eqs. ([3]) or ([4])). 
This is equivalent to our trec,di ff > iff that we discussed 
in §5.3. 

Note that the above suggestion to explain Larson 
(1981) laws is not unique. For instance, there are ways to 
explain Eq. appealing to the properties of MHD tur- 
bulence. Indeed, for the GS95 turbulence, the turbulent 
velocity v scales with the fourth root of the parallel size 

of the eddies, i.e. v ~ Z^^^ and if, as we discussed in §2, 

the parallel scale of turbulent motions is associated with 
the cloud size, Eq. ([25]) follows trivially. Another way to 
reproduce the same law (i.e. Eq. ([25]) ) is to assume steep 
supersonic turbulence with E ~ in clouds, which also 
would result in a similar scalinj^. However, using the 
latter two approaches it is not easy to reproduce the sec- 
ond relation given by Eq. 1261 which decreases the value 
of such derivations. An interesting suggestion was also 
made recently by Kritsuk & Norman (2012). One cannot 
exclude that, in reality, different mechanisms may act to- 
gether providing the Larson (1981) scalings for a variety 
of clouds. 

At the same time, the results in Robertson & Goldrcich 
(2012) indicate that in order to explain the Larson (1981) 
laws appealing to clouds collapse, does require that the 
infall velocities be comparable with turbulent ones. Ob- 
servations (see Tafalla et al. 1998) indicate that the cloud 
collapse is usually slower that that. 

6.2. Expectations for the reconnection diffusion process 

^ The proportionality of the two quantities is valid for super- 
Alfvenic turbulence as well (see Beresnyak 2012). 

We may notice the dependences of the scales on densities ob- 
tained in §5.2 are liower > C\/p~^ and lupper < C2/ p~^J\ where 
Ci and C2 are coefficients. The Larson law given by Eq. I I26II corre- 
sponds to I ~ p~^, and it fits within boundaries determined by the 
upper and lower constraints on the domain of efficient reconnection 
diffusion. 

Steep spectra of turbulence correspond to the observations of 
cold HI in Chepurnov et al. (2010). 



The calculations presented in Santos-Lima et al. 
(2010) testify that reconnection diffusion can naturally 
explain the poor correlation of magnetic field and density 
reported for the diffuse media in Troland & Heiles (1986). 
This result is rather obvious from the physical picture of 
reconnection diffusion that we presented above. Indeed, 
eddies with different magnetization depicted in Figure [3] 
are expected to exchange partially ionized gas irrespec- 
tively of the degree of gas ionization. Therefore the high 
ionization of media that was dealt with by Troland & 
Heiles (1986) should not present a problem for decorre- 
lating densities and magnetic fields. Figure [5] testifies 
that the decorrelation extends to higher densities and, 
moreover, H I and lower density molecular clouds tend 
to have the same magnetization (equal to that in the 
large-scale interstellar medium) that is independent of 
their density. This corresponds to the expectation of the 
reconnection diffusion procesfQ- 

In Lazarian & Vishniac (2009) it was suggested that 
reconnection diffusion can remove magnetic fields from 
protostellar accretion disks to explain observations dis- 
cussed in Shu et al. (2005). Detailed MHD calculations 
in Santos-Lima et al. (2012) showed that the process 
can resolve so called "magnetic breaking catastrophe" of 
the disks. Another problem where reconnection diffusion 
may be valuable is the so called "magnetic flux problem" 
(see Galli et al. 2006 and references therein). For exam- 
ple, T-Tauri stars have magnetic fields 2 x 10"^ Gauss 
(see Johns-KruU 2007) or $/Af sa 3x10"* gauss cm^ g-i, 
which is a million times smaller that the fiux to mass 
ratio estimated for a one solar mass clump in a cloud 
of density 10^ cm~^ (see Draine 2011). While further 
detailed quantitative modeling of the problem would be 
valuable, the numerical results obtained with reconnec- 
tion diffusion are encouraging. For instance, the removal 
of magnetic flux via reconnection diffusion was reported 
both at the stage of turbulent cloud formation (Santos- 
Lima et al. 2010, Leao et al. 2012) and at the stage of 
the accretion disk evolution in the presence of turbulence 
(Santos-Lima et al. 2012). 

It is encouraging that the reconnection diffusion con- 
cept provides new ways of approaching challenging 
physics of star formation in extreme environments. For 
instance, galaxies emitting more than 10^^ solar lumi- 
nosities in the far-infrared are called ultra-luminous in- 
frared galaxies or ULIRGs. The physical conditions in 
such galaxies are extreme with a very high density of 
cosmic rays (see Papadopoulos et al. 2011). Ambipolar 
diffusion is expected to be suppressed due to cosmic ray 
ionization. At the same time these environments have 
the highest star formation rates, which is suggestive of 
a process that removes magnetic fields independent of 
the level of ionization. Reconnection diffusion is such a 
process. 

There are more cases when the predictions based on 
reconnection diffusion provide solutions to long standing 

An alternative explanation of the Troland & Heiles (1986) 
results was provided in Passot & Vazquez-Semadeni (2003) and it 
is based on the notion that MHD turbulence is a superposition of 
turbulent cascades of Alfven, slow and fast modes (Cho & Lazarian 
2002, 2003) and therefore the enhancements of magnetic field and 
densities may correlate or anti-correlate within different modes. 
Unlike reconnection diffusion, this cannot explain the results in 
Figure [S] 
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problems. For instance, the observed star formation rate 
is about the same in galaxies with low metallicities as 
in galaxies with high metallicities (Elmegreen & Scalo 
2004). This is hard to understand if magnetic field loss 
is governed by ambipolar diffusion, which is supposed 
to be much faster in low-Z galaxies with high metallic- 
ities. As a result, contrary to observations, ambipolar 
diffusion theory predicts that the Initial Mass Function 
(IMF) and therefore the star formation rate would be 
shifted in \ow-Z galaxies. As we discussed earlier, recon- 
nection diffusion is independent of metallicity, and thus 
provides a solution consistent with observations. 

In addition, an important property of reconnection dif- 
fusion is that, unlike ambipolar diffusion, it not only pro- 
vides the removal of the magnetic field, but also induces 
its turbulent mixing, which tends to make the distribu- 
tion of magnetic field uniform. Turbulent mixing helps 
in keeping the diffuse media in a magnetized subcritical 
state. In this situation, the external pressure is impor- 
tant for initiating collapse, which well corresponds to the 
observations of numerous small dark clouds not form- 
ing stars in the inter-arm regions of galaxies (Elmegreen 
2011). 

Finally, the reconnection diffusion concept implies that 
in realistic turbulent media there is no characteristic 
density for the collapse to be initiated. Therefore any 
cloud with the appropriate virial parameter (see McKee 
& Zweibel 1992) can form stars. The difference between 
different clouds stems from the density controlling the 
timescale of the collapse and the level of steering pro- 
vided by turbulence or temperature of the media. Di- 
rectly, the requirement of clouds to be molecular is not 
present for the reconnection diffusion to induce star for- 
mation. However, molecular clouds have lower tempera- 
ture. 

6.3. Mean magnetic field and reconnection diffusion 

The reconnection diffusion process, as any diffusion 
process, is expected to relax the sharp local changes of 
magnetic field direction. However, as reconnection diffu- 
sion is associated with turbulence, the latter introduces 
variations of magnetic field direction. 

In terms of reconnection diffusion some observations 
are definitive, some are suggestive. For instance, we be- 
lieve that reconnection diffusion can explain the align- 
ment of magnetic field of molecular clouds cores with the 
magnetic field of the spiral arms (see the example of M33 
measurements in Li & Henning 2011). Indeed, it seems 
natural that if molecular clouds are self-gravitating, the 
gravitational contraction associated with the cloud for- 
mation is expected to distort the local direction of mag- 
netic field even in the case of subAlfvenic turbulence 
of the ordered galactic magnetic fields (see Eq. ([T])). 
Similarly, the distortions are expected to arise from the 
molecular clouds sharing its angular momentum with the 
ambient media during their contraction. Whatever the 
corresponding Alfvenic Mach number Ma the reconnec- 
tion diffusion can smooth and relax the structure of the 
otherwise distorted mean magnetic field. We have the 
extensive evidence of the molecular clouds being turbu- 
lent (see McKee & Ostriker 2007 and ref. therein) and 
therefore the process of reconnection diffusion there is 
inevitable. This process will tend to move magnetic in 
molecular clouds correlated with the mean magnetic field 



removing the distortions induced by the formation pro- 
cesf[3- An alternative interpretation of the results of Li 
& Henning (2011) is that the interstellar matter is col- 
lected along magnetic field lines weakly distorted by sub- 
Alfvenic turbulence. While the quantitative analysis of 
the latter suggestion is beyond the scope of the present 
work, we would like to stress that the magnetic recon- 
nection is arguably an essential part of subAlfvenic mag- 
netic turbulence (see LV99, Eyink et al. 2011). There- 
fore numerical simulations of the accumulation of matter 
along magnetic field lines within subAlfvenic turbulence 
do intrinsically include the magnetic field smoothing by 
reconnection diffusion. While we accept that without a 
detailed study the observations by Li & Henning (2011) 
do not provide a conclusive test of reconnection diffusion, 
we feel that reconnection diffusion provides a viable ex- 
planation consistent with the observations. 

Another possible explanation of molecular clouds hav- 
ing magnetic field correlated with the mean galactic mag- 
netic field is related to the accumulation of matter along 
magnetic field lines. With our paper being focused on 
the comparison of the predictions of reconnection versus 
ambipolar diffusion we leave a detailed comparison of 
the reconnection diffusion with the aforementioned non- 
diffusive scenario of matter collection for further publi- 
cations. We may mention parenthetically that reconnec- 
tion diffusion is expected to affect significantly such a 
collection (see Appendix D). 

7. DISCUSSION 
7.1. Reconnection dijjusion: simulations and estimates 

This paper was motivated by the observational results 
reviewed by Crutcher (2012), who showed that several 
important characteristics were not compatible with the 
ambipolar diffusion theory. We have shown that recon- 
nection diffusion can produce agreement with the obser- 
vations. Reconnection diffusion is based on the numeri- 
cally tested LV99 model of fast reconnection in turbulent 
media as well as more recent insights into the violation 
of the frozen-in condition in turbulent magnetized fluids 
(Eyink 2011, ELVll). It is also appeals to the direct 
3D MHD simulations of the magnetic field diffusion in 
Santos-Lima et al. (2010, 2012) that showed consistency 
with the statistical description of the magnetic diffusion 
in turbulent fluids (see Lazarian 2006). Therefore in the 
present paper we use this statistical description directly. 
This allows us to get insight into the evolution of mag- 
netic fields at the core and the envelope. 

The approach to studying reconnection diffusion that 
is based on 3D numerical simulation has its own limi- 
tations. In Santos-Lima et al. (2010) we used an exter- 
nal gravitational potential. In more recent simulations in 

^•^ In view of the discussion above an additional explanation is 
required. Why would not reconnection substantially distort the 
original direction of the magnetic field? Indeed the visualization 
of simulations of the reconnection in the presence of turbulence 
(see Lazarian, Eyink & Vishniac 2012) show substantial changes of 
the magnetic field direction introduced by magnetic reconnection. 
However, these changes are in the component where magnetic field 
reverses its direction. In the case of reconnection diffusion we are 
dealing with the magnetic fields that may be originally at a small 
angle to each other (see Figures [2][3ll . It is evident from topolog- 
ical considerations that reconnection of such fields cannot create 
a substantial change of the observed magnetic field direction (see 
also the discussion in §2). 
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Leao et al. (2012) self-gravity is used. However, using an 
isothermal code to make the turbulent cloud more well- 
defined in the absence of the confining pressure of the 
ambient gas at a different temperature, we imposed an 
additional gravitational potential. We are experimenting 
with the ideas of introducing the ambient pressure with- 
out compromising the resolution and the simplicity of the 
simulation interpretation, but this is work in progress. 
In this situation, it is advantageous to get estimates and 
analyze idealized situations with clear physical meaning 
the way we do in this paper. Naturally, our simplified 
calculations do not seek the quantitative agreement with 
the observations, but only illustrate the qualitative be- 
havior. Indeed, within the adopted toy model we do not 
evolve density together with magnetic field, but consider 
the diffusion of magnetic field out of an idealized cloud 
of constant densitjo. 

We can mention that it is possible to argue that the 
process of reconnection diffusion is present in some of 
high resolution numerical simulations. However, with- 
out clear identification of the role of turbulence in fast 
reconnection, one may not be sure when the results are 
due to physically motivated reconnection diffusion and 
when they are the consequence of the bogus effects of 
numerical diffusion. For instance, Crutcher, Hakobian 
& Troland (2009) refer to the simulations in Luntilla et 
al. (2009) that produce, in agreement with observations, 
higher magnetization of the cloud cores. If these cores 
are of the size of several grid units across, numerical ef- 
fects rather than reconnection diffusion may be dominant 
and turbulence is suppressed at these scales. 

Our results show that in the presence of reconnection 
diffusion it is natural to expect that the magnetic field 
diffuses faster from the area of the envelope and the mag- 
netic field strength is larger at the cloud core. This con- 
clusion does not really depend on the particular model 
of turbulence and its distribution in the molecular cloud 
as far as those are constrained by existing observations. 

We should also stress that the 3D numerical simu- 
lations like those in Santos-Lima et al. (2010, 2012) 
by themselves do not provide the description of recon- 
nection diffusion in realistic astrophysical environments. 
The interpretation of the results requires the proper un- 
derstanding of scaling of magnetic reconnection with 
the dimensionless combination called the Lunquist num- 
ber S = {Lcur.shVA/'n)i where Lcur.sh. is the extent of 
the relevant current sheet, is Ohmic diffusivity. The 
Lundquist numbers in molecular clouds and in the corre- 
sponding simulations differ by a factor larger than 10^. 
In this situation one can establish the correspondence 
between the numerical simulations and astrophysical re- 
ality only if the reconnection does not depend on S. The 
independence from S of magnetic reconnection is the con- 
clusion of LV99 model. This model has been tested in 
Kowal et al. (2009, 2012) via a set of dedicated numeri- 

A possible initial set up with the cloud of gas density ~ 60 
cm~^, radius 2.5 pc and the initial magnetic field of ~ 6 X 10~^G 
may be advantageous. Such a cloud would be subcritical, which 
agrees with the notion of diffuse H I clouds being subcritical. Then 
the reconnection diffusion would remove flux such that the cloud 
could collapse, with density building up to 10'' cm~^ in the r = 0.2 
pc core and up to 10^ cm~^ in the r = 1 pc envelope, which are the 
mean parameters of the clouds observed by Crutcher, Hakobian & 
Troland (2009). This is what we will try in future. 



cal simulations that confirmed the scaling predictions in 
LV99. This work exemplifies the advantages of the syn- 
ergy of scaling arguments with numerical simulations as 
opposed to the "brute force numerical approach" , which 
may not be productive while dealing with turbulence. 

7.2. Predictions of the reconnection and ambipolar 

diffusion models 

Within the concept of ambipolar diffusion the explana- 
tion of observational results reviewed by Crutcher (2012) 
is extremely difficult. Our point is that ambipolar diffu- 
sion is not the only process that can be responsible for 
the removal of magnetic flux from the media. We show 
that reconnection diffusion can reduce the magnetic field 
strength in the envelope with respect to the core. On a 
more fundamental level, reconnection diffusion changes 
the mass to flux ratio allowing magnetic fleld to diffuse 
away from the center of the gravitational potential. 

The available observations provide a rather mild con- 
straint, namely, that mass/flux does not increase as fast 
in the core as ambipolar diffusion predicts. In many 
ways, ambipolar diffusion is a very special type of diffu- 
sion, the efficiency of which drops towards more ionized 
outer regions of the envelope inducing more efficient flux 
loss from cores. 

Typically the observed column densities through the 
cores studied by Crutcher et al. (2009) are about twice 
those through the envelope regions. If mass/flux were 
constant, then the fields in the envelopes would be 1/2 
those in the cores. So just saying that the fields in the 
envelopes are less than those in the cores does not require 
reconnection diffusion, or any kind of diffusion. Recon- 
nection diffusion is required if the envelope fields are even 
weaker than half of the core fields; this is what is ob- 
served. Ambipolar diffusion would require that envelope 
fields are stronger than half of the core fields. 

7.3. Reconnection diffusion and ambipolar drift in 

turbulent media 

Within this study we do not appeal to ambipolar diffu- 
sion, which is acceptable when the reconnection diffusion 
provides larger diffusivity for magnetic fields. One may 
argue that this is a generic situation in the presence of 
turbulence. 

For instance, Heitsch et al. (2004, henceforth HX04) 
performed 2.5D simulations of turbulence with two-fluid 
code and examined the decorrelation of neutrals and 
magnetic field that was taking place as they were driving 
the turbulence. The study reported an enhancement of 
the ambipolar diffusion rate compared to the ambipolar 
diffusion acting in a laminar fluid. HX04 correctly asso- 
ciated the enhancement with turbulence creating density 
gradients that are being dissolved by ambipolar diffusion 
(see also Zweibel 2002). While in 2.5D simulations of 
HX04 the numerical set-up precluded reconnection from 
happening (as magnetic field was perpendicular to the 
mixing plane and magnetic fields never crossed each other 
at an angle) , the authors reported an enhanced rate that 
is equal to the turbulent diffusion rate LVl, which is 
the result expected as a limiting case of the reconnection 
diffusion prediction for the special set-up studied. 

While we agree with HX04 in terms of the importance 
of turbulence, we suggest that it is misleading to refer 
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to the process as "turbulent ambipolar diffusion". We 
believe that ambipolar diffusion does not play a role in 
the enhancement measured and the observed effect is en- 
tirely due to turbulenc43- We note that, in the presence 
of turbulence, the independence of the gravitational col- 
lapse from the ambipolar diffusion rate was reported in 
numerical simulations by Balsara, Crutcher & Pouquet 
(2001). 

Ambipolar diffusion may control flux removal from the 
laminar cores. In the presence of turbulence, it is recon- 
nection diffusion that dominates. The results of HX04 
support this notion, as the authors in their set up re- 
port the transition to turbulent diffusivity in the pres- 
ence of turbulence. The set-up in HX04 precludes the 
magnetic fields from reconnection as field lines perform 
mixing motions being kept absolutely parallel. In any 
realistic turbulent 3D flow magnetic reconnection will be 
essential. 

7.4. Reconnection diffusion and hyper-restivity concept 

Another process of diffusive nature is related to "hyper- 
resistivity" or enhanced physical resistivity of turbulent 
fluids. To explain fast removal of magnetic field from 
accretion disks Shu et al. (2006) appealed to the hy- 
perrestivity concept (Strauss 1986, 1988, Bhattacharjee 
& Hameiri 1986, Hameiri & Bhattacharjee, Diamond & 
Malkov 2003). The studies introducing hyper- resistivity 
attempt to derive the effective resistivity of the turbulent 
media in the context of the mean-field resistive MHD. 
Using magnetic helicity conservation the authors derived 
the electric field. Then, integrating by part, they ob- 
tained a term which could be identified with effective 
resistivity proportional to the magnetic helicity current. 
There are several problems with this derivation. In par- 
ticular, the most serious is the assumption that the he- 
licity of magnetic field and the small scale turbulent 
fields are separately conserved, which erroneously dis- 
regard the magnetic helicity fluxes through open bound- 
aries that is essential for fast stationary reconnection (see 
ELVll). In more general terms, hyper-resistivity idea is 
an incarnation of the mean-field approach to producing 
fast reconnection. As explained in ELVll, the problem 
of such approaches is that the lines of the actual astro- 
physical magnetic field should reconnect, not the lines of 
the mean field. Therefore the correct approach to fast 
reconnection should be independent of spatial and time 
averaging. 

All in all, we believe that the concept of hyper- 
resistivity should not be applied to astrophysical envi- 
ronments. 

7.5. Reconnection diffusion and "magnetic turbulent 
diffusivity" concept 

A similar process takes place in the case of molecular diffu- 
sivity in turbulent hydrodynamic flows. The result for the latter 
flows is well known: in the turbulent regime, molecular diffusivity 
is irrelevant for the turbulent transport. Indeed, in the case of high 
microscopic diffusivity, turbulence provides mixing down to a scale 
at which the microscopic diffusivity both suppresses the cascade 
and ensures efficient diffusivity of the contaminant. In the case of 
low microscopic diffusivity, turbulent mixing happens down to a 
scale I2 <SC li, which ensures that even low microscopic diffusivity 
is sufficient to provide efficient diffusion. In both cases the total 
effective diffusivity of the contaminant is given by the product of 
the turbulent injection scale and the turbulent velocity. 



One may claim that the reconnection diffusion concept 
extends the concept of hydrodynamic turbulent diffusion 
to magnetized fluids. The physics of it is very different 
from the "magnetic turbulent diffusivity" idea discussed 
within the theories of the kinematic dynamo (see Parker 
1979). Reconnection diffusion, unlike "magnetic turbu- 
lent diffusivity" , deals with dynamically important mag- 
netic fields, e.g. with subAlfvenic turbulencc[3. Thus in 
the process of reconnection diffusion, magnetic fields are 
not passively mixed and magnetic reconnection plays a 
vital role for the process. 

One may show that the domain of applicability of 
the "magnetic turbulent diffusivity" concept is extremely 
limited. Numerical calculations (Cho et al. 2010, Beres- 
nyak 2012) testify that around 5% of turbulent energy 
fiow is transferred to magnetic energy and that this value 
does not depend on the initial magnetization of the fluid. 
Therefore, even if the turbulent flow initially is not mag- 
netized, it is expected to develop dynamically impor- 
tant magnetization on the time scale of 10-20 large eddy 
turnover times. As a consequence, we expect that star 
formation will happen in magnetized media even for the 
first generation of starsFl 

7.6. Reconnection diffusion and numerical effects 

One might argue that reconnection diffusion is auto- 
matically a part and parcel of 3D numerical simulations 
of star formation. As we argue below this is only par- 
tially true and a proper understanding of reconnection 
diffusion is required for the interpretation the numerical 
simulations. 

The LV99 model predicts that the reconnection rates in 
turbulent fluids are independent of the local physics, but 
are determined by the turbulent motions. This is great 
news for numerical simulations: parasitic numerical ef- 
fects that induce poorly controlled small scale diffusivity 
of magnetic field lines are not important on the scales 
of turbulent motions. In other words, the low resolu- 
tion numerics may provide an adequate representation 
of high Lundquist number astrophysical turbulence as 
far as the reconnection diffusion is concerned. On the 
contrary, if the structures studied in numerical simula- 
tions (e.g. cores, filaments, shells) lose turbulence due 
to numerical diffusivity effects, we predict that the diffu- 
sion of magnetic flux in those simulated structures differs 
signiflcantly from the diffusion in the actual interstellar 
structures where turbulence persists. Note, that naive 
convergence studies based on increasing the numerical 
resolution several times cannot notice the problem unless 
the resolution increased to the degree that the aforemen- 
tioned structures become turbulent. 

In terms of numerical studies of star formation it is 
encouraging that it is the largest scales of turbulent mo- 
tions that are important for the reconnection diffusion. 
At the same time, in many cases in numerical simula- 

The concept of "magnetic turbulent diffusivity" assumes that 
magnetic fields can be passively mixed up on all scales, up to the 
Ohmic diffusion scale. At the latter scale Ohmic dissipation is 
trivial and the diffusion is hydrodynamic at all scales. 

Incidentally, the same logic suggests that the proposals of 
bringing magnetic field and kinetic energy to a state of equipar- 
tition in molecular clouds during their formation (see Sur et al. 
2012) may not work. Indeed, molecular clouds are not likely to 
survive 20 crossing times. 
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tions we expect to observe magnetic field evolution influ- 
enced by strongly damped turbulent motions, which do 
not demonstrate power law cascade. The differences of 
reconnection diffusion in such damped regimes and in the 
astrophysically important regime of well developed tur- 
bulence with a power law spectrum should be a subject 
of further studies. 

At the same time, the importance of understanding re- 
connection diffusion is surely not limited by the proper 
interpretation of numerical simulations of star formation. 
Reconnection diffusion provides a theoretical way to pre- 
dict the dynamics of magnetic flux in turbulent fluids. 

7.7. Magnetic star formation in the presence of 
reconnection diffusion 

In spite of the fact that reconnection diffusion relaxes 
the conventional flux freezing condition, the role of mag- 
netic fields for star formation should not be disregarded. 
For instance, magnetic fields solve the angular momen- 
tum problem in the process of star formation (Mestel 
1965). At the same time, reconnection diffusion allows to 
prevent the "magnetic breaking catastrophe" that, con- 
trary to observations, prevents the formation of circum- 
stellar disks (Lazarian & Vishniac 2009, Santos-Lima et 
aL 2012). 

7.8. Comparison with earlier papers 

The concept of removal of magnetic fields from clouds 
via the process of LV99 reconnection was first introduced 
in Lazarian (2005) (see also discussions in Lazarian & 
Vishniac 2009). The numerical calculations of reconnec- 
tion diffusion in application to molecular clouds and ac- 
cretion disks have been performed in Santos-Lima et al. 
(2010, 2012). However, the physical picture of how mag- 
netic fields and matter can decouple in turbulent magne- 
tized plasma was not presented there. 

In comparison, this paper focuses on the physics of 
reconnection diffusion and presents a simple analytical 
model illustrating how the tendencies in magnetic field 
distribution revealed in the observations by Crutcher et 
al. (2009, 2010a) can be accounted for, as well as other 
observations reviewed by Crutcher (2012). In addition, 
the current paper presents the additional evidence of the 
importance of reconnection diffusion in various astro- 
physical environments. While alternative explanations 
may be possible for the cases discussed (e.g, Li et al. 
2011) we feel that the consistency of those with the recon- 
nection diffusion concept provides an additional evidence 
in favor of the reconnection diffusion scenario. This calls 
for the necessity to include the concept of reconnection 
diffusion within the star formation paradigm. 

8. SUMMARY 

This paper presents the physical justifications of the 
concept of reconnection diffusion and an analytical toy 
model as well as simple numerical calculations of the 



effects of the process on the magnetic fiux loss out of 
clouds. The major results of this paper can be briefly 
summarized in the following way: 

• Reconnection in turbulent media happens fast, and 
this enables turbulent diffusion of matter and mag- 
netic field. This process is termed reconnection dif- 
fusion, and it provides a way to change the fiux to 
mass ratio of molecular clouds. 

• The process of reconnection diffusion does not de- 
pend on the ionization of the cloud material; there- 
fore, if the magnetic field is captured within a col- 
lapsing turbulent cloud, the decrease of the mag- 
netic field strength due to reconnection diffusion is 
most prominent at the edges of the cloud. This is 
a robust result that does not depend on the details 
of the turbulence picture. 

• The predictions of reconnection diffusion are con- 
sistent with the measurements of the ratio of the 
magnetization of cores and envelops of molecular 
clouds obtained in Crutcher et al. (2009), with 
the poor density - magnetic field correlation at low 
densities, with the 2/3rds power law correlation at 
high densities reported in Crutcher et al. (2010b), 
and with the large range in magnetic field strengths 
in molecular clouds of the same density. 

• In agreement with observations, reconnection dif- 
fusion predicts the value of critical column density 
at which the increase of magnetic field with den- 
sity is measured. Reconnection diffusion may also 
provide an explanation of the well-known empirical 
Larson relations. 

• Reconnection diffusion can provide a way of ap- 
proaching other problems where ambipolar diffu- 
sion is unlikely to produce viable solutions, e.g. 
the problems (to name a few) of poor correlation of 
density and magnetic field in diffuse media, efficient 
star formation in ultra-luminous infrared galaxies, 
alignment of magnetic field in cloud cores with the 
magnetic field in spiral arms. 

• We do not argue here that reconnection diffusion 
is the only important MHD process that may af- 
fect star formation. If subcritical, self-gravitating 
clouds are formed that have weak turbulence, am- 
bipolar diffusion would operate over longer time 
scales to drive the cores to supercritical states and 
allow stars to form. It is also possible that collec- 
tion of matter along magnetic flux tubes play a role, 
perhaps especially at smaller spatial scales where 
GMCs fragment to form individual self-gravitating 
clouds. It seems likely that turbulence would pre- 
vent the requisite laminar flow, and reconnection 
diffusion would come into play. 



APPENDIX 

A: DERIVATION OF LV99 RECONNECTION RATE 

A picture of two flux tubes of different directions which get into contact in 3D space is a generic framework to 
describe magnetic reconnection. The upper panel of Figure [1] illustrates why reconnection is so slow in the textbook 
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Sweet-Parker model. Indeed, one can write 

A 

Wrec = VA — , (Al) 

meaning that Vrec ^ Va if A ^ L^;. Consider the Ohmic diffusion of magnetic field lines (see ELVll). The mean-square 
vertical distance that a magnetic field line can diffuse by resistivity in time t is 

{y\t)) ^ \t, (A2) 

where A = -qc? /Air is Ohmic diffusivity. The field lines are advected out of the sides of the reconnection layer of length 
Lx at a velocity of order Va- Thus, the time that the lines can spend in the resistive layer is the Alfven crossing time 
Ia = Lx/Va- Thus, field lines can only be merged that are separated by a distance 

A = ~ \/A<I = Lx/yfS, (A3) 

where S is Lundquist number, 

S = LxVa/\. (A4) 
Combining Eqs. (|A1[) and (jA3l) one gets the famous Sweet-Parker reconnection rate, 

Vrec,SP ^Va/^- (A5) 

The LV99 model generalizes the Sweet-Parker one by accounting for the existence of magnetic field line stochasticity 
(Figure [T] [lower panels]). In ELVll the LV99 reconnection rate was rederived appealing to the concept of Richardson 
(1926) diffusion. Richardson diffusion (see Kupiainen et al. 2003) implies the mean squared separation of particles 

{\x^{t) - X2{t)\') ^ et\ (A6) 

where t is time, e is the energy cascading rate and () denote an ensemble averaging. For subAlfvenic turbulence 
e K u jj (VaL), where ul is the injection velocity and L is an injection scale (see LV99). Therefore analogously to Eq. 
(IA3p . one can write 

A « « Lx{Lx/L)'/H'll, (A7) 

where it is assumed that Lx < L. Combining Eqs. (|A1|) and (|A7|) . one recovers the LV99 expression for the rate of 
magnetic reconnection (see also ELVll), 

Vrec,LV99^VAiLx/L)^/^Ml, (A8) 

in the limit of Lx < L. Analogous considerations allow one to recover the LV99 expression for Lx > L, which differs 
from Eq. ^ by the change of the power 1/2 to —1/2. 

B: MICROSCOPIC PICTURE OF RECONNECTION DIFFUSION 

Magnetic field wandering provides a possibility of explaining reconnection diffusion without appealing to magnetic 
eddies. According to LV99, the transverse displacement of magnetic field lines over distance a; is a spatial random 
walk given by the equation d{y'^)/dx ^ L{Vl/Va)'^ , which results in 

(y2) ^ Lx{VL/VAf. (Bl) 

Figure [TOl left, illustrates the spread of magnetic fie ld li nes in the perpendicular direction as magnetic field lines are 
traced by particles moving along them. Using Eq. (|B1|) one can get the RMS separation of the magnetic field lines 
(LV99) 

i.e. is proportional to /jj. This regime identified in LV99 is a Richardson diffusion but in terms of magnetic field 

lineJ*^. The numerical testing of this prediction of LV99 is presented, e.g. in Lazarian et al. (2004). This is the 
case of diffusion in space. Eq. (jB2|) allows one to calculate the distance l\\ at which the magnetic field lines of regions 

separated by l± start overlapping^, i-e. /|| ~ {61\LY/'^{Va/Vl)'^/^ . Thus for sufficiently large l\\ all parts of the volume 
of magnetized plasmas get connected. In other words, the entire volume becomes accessible to particles moving along 
magnetic field lines. Naturally, in this situation the star formation community's customary notion of flux to mass ratio 
loses its original sense. 

In addition, in turbulent plasmas, our turbulent volumes spread due to Richardson diffusion with 6l\ ~ . This 
process is illustrated by Figure IIOI right. The magnetized plasma spreads by subAlfvenic turbulence over a larger 

This regime induces perpendicular "superdifTusion" in terms 19 t -j i n ■ -t i -t- mnl ii, 

„ . I J, I , J. 1 i -1 Incidentally, m terms oi reconnection, Eq. I1B2JI expresses the 

ol cosmic rays and other charged particles streaming along mag- ■ i r iu j-a j * j a ■ t^- c u j- i i- 

.„,,..■' Of o b b thickness ol the o utflo w denoted as A m figure [IJ Substituting 

netic held lines. ^j^j^ mine in Eq. i EH i. one recovers the LV99 reconnection rate, 

i.e. Eq. JlJ. 
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51 I ~ t5 



Case A: in space 



Case B: in time 




Fig. 10. — Microscopic physical picture of reconnection diffusion. Magnetized plasma from two regions is spread by turbulence and mixed 
up over Ai„t- Left panel: Description of the process in terms of field wandering in space. Right panel. Description of the magnetic field 
line spread with time. 

volume, while motions of plasma along magnetic field lines are ignored. The two initially disconnected volumes get 
overlapped, and the magnetic fields and plasmas of two regions get mixed up. Again, one easily can see that the 
process that we described in terms of magnetic field lines is similar to the one we described in terms of reconnected 
magnetic field fiux tubes in Figure [3l This is expected, as magnetic reconnection is an intrinsic part of the picture of 
MHD turbulence that governs the dynamics of mag netic field lines (LV99j3 

C: TRACING MAGNETIC FIELD LINES WITH IONS 

Star formation theory was developed in the assumption that ions trace magnetic fields. Thus two ions on the same 
magnetic field line are expected to stay on the same line unless their motion is perturbed by collisions. This notion 
needs to be modified in the presence of turbulence. 

The study in LV99 revealed the Richardson-type diffusion of magnetic field lines. Figure [11] illustrates the loss of 
the Laplacian determinism for magnetic field lines (see also Eq. [TO|l . In analogy with the illustrative example above, 
the final line spread l± does not depend on the initial separation of the field lines. This is a remarkable effect that 
provides a microscopic picture of reconnection diffusion based on the description of magnetic field lines rather than on 
the reconnection of well-organized fiux tubes. 

We shall consider tracing magnetic field lines in the realistic turbulence situation with the dissipation scale lmin,±i 
where, as everywhere in the paper, _L denotes the scale perpendicular to the local magnetic field (see Figure[Tl]). While 
a detailed discussion of the Richardson type diffusion in the presence of magnetic field is given in ELVll, we feel that 
our introduction of the finite dissipation scale lmin,± helps to avoid some of the paradoxes discussed in ELVll. It is 
also important that for the problems of star formation lmin,± may be much larger than ion Larmor radius. In other 
ways we follow the logic presented in ELVll. First of all, resistivity, whatever its nature, introduces stochastic forcing 
in the description of magnetic field line dynamics. Indeed, the induction equation with the resistive term rjAB signifies 
stochasticity associated with Ohmic diffusion. Therefore the definition of the magnetic field line on scales affected by 
resistivity is not deterministic. In addition, it was stressed in ELVll that the magnetic field line motion is a concept 
defined by convention and not testable experimentally. This point was discussed in the literature (see Newcomb 1958, 
Vasylianas 1972, Alfven 1976), but sometimes forgotten. Thus magnetic field lines may be tagged by ions that start at 
the same field line (see Figure [TT|) . In the case of smooth laminar magnetic field and ideal MHD equations the motions 
of ions will reveal magnetic field lines and two ions on the same field line will always remain on the same line. The 
situation is radically different in the presence of turbulence. 

In our thought experiment we shall trace ions moving with the same velocity and separated perpendicular to magnetic 
field by a distance of a Larmor radius po. As we discuss further, if the separation is less po than that one can appeal 
to plasma effects to increase the separation to po- Let us assume that the minimal scale of turbulence lmin,± > Po- In 
this situation the dynamics of ions can be approximated by the dynamics of charged particles in toy model of "a single 
scale MHD turbulence" discussed in Rechester & Rosenbluth (1978). Indeed, the turbulent motions at the critically 
damped scale lmin,± are dominant for shearing and steering matter and magnetic field on the smaller scales. The 
Rechester & Rosenbluth (1978) theory predicts the Lypunov growth of the perpendicular separation of ions, i.e. the 
separation gets po exp(^/Zmi„ y, where I is the distance traveled by ions and lmin,\\ is the scale parallel scale of the 
critically damped eddies with the perpendicular scale lmin.± (see also Narayan & Medvedev 2001, Lazarian 2006). 



Thus to get separated by the distance l„ 



the ions should travel the so-called Rechester-Rosenbluth distance 



^RR 



^min,\\^'^{^niin,± / Pq\ 



(CI) 



which is at most a dozen times larger than the microscopic scale lmin.\\- As soon as ions get separated over the 
distance of lmin.±, they get into different eddies and the process of Richardson diffusion starts. The separation of ions 



One can also claim that spontaneous stochasticity of mag- 
netic fields in turbulent fluids is an underlying process that gov- 



erns magnetic reconnection (ELVll). As we mentioned earlier, fast 
reconnection makes MHD turbulence theory self-consistent. 
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Fig. 11. — Left Panel: The extent of th e di fEusion of magnetic field hnes is independent of their initial position in agreement with the 
concept of Richardson diffusion (see Eq. mOt *). Right Panel: Upper plot Ions tracing the same magnetic field line. The diffusion and 
decorrelation arises from plasma or Ohmic effects as well Rechester-Rosenbluth effect. Middle plot: Ions separated by scales much larger 
than the ion Larmor radius are further separated by the Rechester-Rosenbluth effect. Lower plot: At scales larger than the turbulence 
damping scale the Richardson diffusion takes over resulting in explosive separation of field lines. 

increase fast with their initial position quickly forgotten (see the solution given by Eq. (fTO)) and the detailed discussion 
in ELVll). Thus after a relatively short period during which ions move in a correlated manner remembering their 
original position, a stochastic regime when the initial vicinity of the ions is completely forgotten takes over. As we 
used the ions as tracers of magnetic field, we can talk about the stochasticity of magnetic field lines (as traced by 
ions). 

D: RECONNECTION DIFFUSION AND COLLECTING OF MATTER ALONG MAGNETIC FIELD LINES 

Another concept that was discussed in the literature is the one based on changing the flux to mass ratio by collecting 
matter along magnetic field lines. Similar to reconnection diffusion, this process can happen within the MHD approx- 
imation and it does not depend on the degree of media ionization (see Vazquez-Semadeni et al. 2011). Incidentally, 
as one can easily observe in Figure 2, the motion along field lines is also a part of the reconnection diffusion concept. 
The difference between the two aforementioned pictures arises, however, from the fact that within the "motion along 
magnetic field lines" concept, flux freezing holds and the fluid moves along field lines that preserve their identity 
through the process. On the contrary, the reconnection diffusion concept is based on fast reconnection arising from 
the LV99 process or, equivalently, on the violation of flux freezing in turbulent flows that arises from fundamental 
properties of Lagrangian dynamics of these flows (see ELVll). 

It is beyond doubt that the collection of matter along magnetic fleld lines is an important process for laminar flows. 
Reconnection diffusion is not applicable for those as it requires turbulence for enabling both fast reconnection and 
diffusion of magnetic fields. However, as we mentioned earlier, observations testify that turbulence is ubiquitous in 
the ISM. The fact that reconnection diffusion does act to smooth magnetic field lines twisted by turbulence does help 
promote motion of matter along field lines in turbulent media. However, it is necessary to investigate the relative 
importance of the concept of moving matter along turbulent magnetic fleld lines adopting the parameters of turbulence 
consistent with observations and appealing to the scalings arising from the widely accepted GS95 model of magnetic 
turbulenceEH- 

While the motion of fluid along laminar magnetic flelds is well defined, the situation is a bit more complicated 
in turbulent fluids. Thinking about the motion of fluids along turbulent magnetic flelds one may visualize both the 
motion along the mean magnetic fleld and the motion along the actual turbulent field. Below we discuss both cases 
for the model cloud that we discussed in Lazarian et al. 2012, §4. 

Consider first to what extent the actual fluid dynamics can be approximated by the motion along the mean field. 
Assume that the material in the envelope is being collected along a cylinder of denv = 1 pc oriented along the mean 



The latter, while not being exact, is so far the best model 
of MHD turbulence that we are aware of (see numerical testing in 
Beresnyak & Lazarian 2010, Beresnyak 2011, 2012). The appli- 



cability of the GS95 scalings to the compressible turbulence was 
demonstrated via numerical simulations (Cho & Lazarian 2002, 
2003, Kowal & Lazarian 2010). 
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magnetic field. Taking the typical density of interstellar medium of tijsm = 20 cm ^, one can estimate the length of 
the flux tube using the mass conservation relation 

denvelope'^envelope ~ ^tube.env^l S AI : (-^-^) 

which for the adopted value rienveiope — 10'^ cm~'^ provides Dtube,env = 50 pc. If interstellar densities are smaller, the 
collection distance Dtube,env only increases. Assuming that the injection scale of interstellar turbulence is L = 100 pc, 
one has to accept that the magnetic field lines undergo Richardson superdifFusion that obeys Eq. (6) in Lazarian et 
al. (2012) (see also ELV12 and ref. therein, Lazarian et al. (2004) for the numerical verification). Ignoring the initial 
separation of the field lines Iq ~ d^nv one can estimate the separation as 

llff « iVl/L)t^ (D2) 

where we took into account that the cascading rate a in MHD turbulence is Vj^/L if turbulence is injected at the 
velocity Vl = 10 km/s, which is larger than the Alfven velocity in the diffuse gas. 

In the scenario of gas being collected along magnetic field lines, the velocities relevant for creating molecular clouds 
are due to turbulence. The maximal turbulent velocity is limited to the injection turbulent velocitjo that we denoted 
Vl- Therefore, the minimal time for the accumulation of the matter for the envelope is 

Combining Eqs. (|D2I) and (jD3|) one gets 

which for the adopted injection scale L = 100 pc and estimated value of Dtube.env provides the separation of the field 
lines of Idiff ~ 30 pc. The fraction of material that actually moves along the the cylinder oriented along the mean 
magnetic field in our case is ~ d^envl^'diff ~ 10^"*. 

In the case of the core one can assume that the mean magnetic field is constricted at the place of the core in 
accordance to magnetic flux conservation: 

B initial d^j^^ i^^j.^ — Bfinald^^^^ (D5) 

and we will use dini.core as the diameter of the flux tube along which the collection of matter parallel to magnetic 
fleld is considered. Adopting the value of Biniuai of 10 /LtG one gets dini^core ~ 0.3 pc. The same arguments that we 
used for the envelope provide for the core the collection scale of 150 pc, which is larger than the assumed scale of the 
turbulence injection. As a result, spreading of the magnetic fleld lines over this scale exceeds the injection scale L and, 
compared to the envelope, an even smaller fraction of matter can be considered to move along mean magnetic fieldp^ 

As the previous discussion shows that the collection of matter along the mean magnetic field may not be tenable for 
the parameters of turbulence that we adopted, one may wonder whether it may be good to talk about collecting matter 
along the stochastic magnetic field lines but assuming that these field lines are static and the diffusion of magnetic 
field during the collection time is negligible. Consider the envelope case, as it presents a less extreme example of 
matter accumulation. It is easy to see that the issue of "static" is violated, as during the time of the collection, 
the magnetic field and fluid in the core performs {ncioud/ninit){denveiope/ L)^/'^ turnovers which is about 17 for the 
numbers adopted. The diffusion of magnetic field from the core over the collection time obeys the Richardson la w, i.e . 
(y2) « V^t^ /L, which for the collection time tmin provides the estimate that coincides with one given by Eq. (|D2p . 
Therefore, the spread of magnetic field lines induced by the diffusion in turbulent media on the timescale of t„iin is 
very large compared to the scale of the envelope and we have to accept that magnetic field diffusion is very important. 
In other words, it is difficult to justify the collection of the matter along passive and static magnetic field lines. It is 
easy to see that the collection of matter within the core appealing to driving matter along turbulent magnetic field 
lines and ignoring the magnetic field line spread is even more questionable than for the case of the envelope. 

We should admit that our arguments above do not account for self-gravity, which may change the ability of magnetic 
field lines to spread at the core, for instance. However, self-gravity cannot be important over the entire extent of flux 
tubes which is ~ 50 pc for the envelope and ~ 150 pc for the core. Therefore we believe that, while matter is moving 
along magnetic fleld lines, the effects of reconnection diffusion cannot be ignored. In fact, due to the pull of gravity 
we expect that matter in the presence of reconnection diffusion will be collected from all sides, rather than exclusively 
from the local direction of the magnetic fleld. Movement of matter along local fleld lines will of course be stronger 
than along other directions. In other words, the motion of matter is three dimensional within turbulent fluids and 
the notion of a one dimensional approximation assumed within the concept of the "matter collection along magnetic 
fleld lines" is of limited applicability if the present day understanding of MHD turbulence is correct. This, in fact, 
corresponds well to the notion of violation of magnetic flux freezing in turbulent fluids (see ELVll). 



In reality only a fraction of 20% of turbulent energy is in 



^■^ Numerical simulations show some correlation between mag- 



the potential motions that can be associated with compressions i- c u j i -i ■ i ai - ■ i -i i -t i 

. ^ n I 1 1- / ^1 T ■ rr^r^r^nw T uctic faeld and velocity m transAlvemc runs (see Kritsuk et al. 

along magnetic held lines (see Cho & Lazarian (2002 ). In the nnn\ tiu- u-j-if iu j i-r c 

.= ° . . ., r .1 \ " . , 2012). ihis by itself cannot be used as an argument in favor of 

text for obtaining the upper limit tor the compression, we consider n i- ■ i j i i ■ i ? n j- i 

,, ,. Til ui 1 i 1 i ■ ■ i- collecting matter in clouds along large msterstellar distances, as 

the limiting case of the compressible velocity equal to the m ection i- c u i- j ■ j j- 

J . " ^ J 1 J magnetic held lines are wandering and spreading. 
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